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Abstract 

The present paper deals with mathematical models of heat and moisture 
transport in layered building envelopes. The study of such processes gener- 
ates a system of two doubly nonlinear evolution partial differential equations 
with appropriate initial and boundary conditions. The existence of the strong 
solution in two dimensions on a (short) time interval is proven. The proof 
rests on regularity results for elliptic transmission problem for isotropic com- 
posite materials. 
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1. Introduction 

Building envelopes, which act as barriers between the indoor and outdoor 
environments, present a crucial component responsible for the building's per- 
formance over the whole service life. In this regard, an important requirement 
to achieve an energy-efficient design is the assessment of the heat and mois- 
ture behavior of the component when exposed to natural climatic conditions. 
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This task can hardly be accomphshed by purely experimental means, mainly 
due to the long-term character of environmental variations and the trans- 
port processes involved. Therefore, a considerable research effort has been 
devoted to the development of predictive models for coupled heat and mass 
transfer in building materials, see e.g. js], lof for historical overviews. 

The major challenge in predicting the transport phenomena in building 
components lies in their complex porous microstructure, resulting in an intri- 
cate mechanism of moisture absorption from surrounding environment. Here, 
the dominant physical processes involve adsorption forces, attracting vapor 
phase molecules to solid parts of the porous system, and capillary conden- 
sation in pores. This needs to be complemented with non-linear dependence 
of thermal conduction on temperature and water content. As a result, engi- 
neering models of simultaneous heat and moisture transfer are posed in the 
form of strongly coupled parabolic system with highly non-linear coefficients. 
Discretization of these equations, typically based on finite volume or finite 
element methods, then provides the basis for numerous simulation tools used 



in engineering practice, see e.g. for a recent survey. However, to our best 
knowledge, the qualitative properties of the resulting systems remain largely 
unexplored. 

The mathematical models of transport processes in porous composite 
media consist of the balance equations, governing the conservation of mass 
(moisture) and thermal energy, supplemented by the appropriate boundary, 
transmission and initial conditions. This system can be written in the form 

^^I^-V.A\{u,,Vu,) = fl{u,) mQ,T, J = 1,2, (1) 
with the nonlinear boundary conditions 

- Al{ui,Vu(,) ■ ne{x) = gj{x,t,ul) onSer, j = l,2, (2) 
the so-called transmission conditions 



(3) 



-Al{ue, Vui) ■ ne{x) = -A'^{Um, Vu^) ■ n^(x) on E^^, 

j = 1,2, and the initial condition 

Ui{x,0) = fi^{x) inQi. (4) 

Here, Q represents a two-dimensional bounded domain with a Lipschitz- 
continuous boundary F = dQ; n = [ni, 77-2) denotes the outer unit normal to 
r. fl consists of M disjoint subdomains fli with boundary dfle, i = I, ■ ■ ■ , M, 
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separated by smooth internal interfaces = dflm H dflg ^ 0. For a fixed 
positive T, we denote by Qit the space-time cyhnder Qit = y< (0,T), 
similarly Ser = (Sl^^nF) x (0, T) and E^^ = Tme x (0, T). Further, in dlD-dl]), 
= (uljuj) represents the unknown fields of state variables and the vector 
t^i = if^lyf^I) describes the initial condition. By Bi, Aj, f^, g^, we denote 
the vectors = iB},B]), Aj = (Af,Af), /, = g, = igj,gj), 

which are functions of primary unknowns ug. Hence, the problem is strongly 
nonlinear. 

The existence of weak solutions to the system ([1]) in homogeneous bounded 
domains [i = 1) subject to mixed boundary conditions with homogeneous 
Neumann boundary conditions has been shown by Alt and Luckhaus in j^. 
They obtained an existence result assuming the operator B in the parabolic 
part to be only (weak) monotone and subgradient. This result has been ex- 
tended in various different directions. Filo and Kacur 0| proved the local 
existence of the weak solution for the system with nonlinear Neumann bound- 
ary conditions and under more general growth conditions on nonlinearities 
in u. These results are not applicable if B does not take the subgradient 
form, which is typical of coupled heat and mass transport models. 

In this context, the only related works we are aware of are due to Vala [26^], 



Li and Sun [20| and Li et al. [21[. Nonetheless, even though [26| admits 
non-symmetry in the parabolic part, it requires unrealistic symmetry in the 
elliptic term. The latter works, studying a model arising from textile in- 
dustry, prove the global existence for one- dimensional problem using the 
Leray-Schauder fixed point theorem. The proofs, however, exploit the spe- 
cific structure of the model and as such are not applicable to our general 
setting. 



In this paper we adapt ideas presented by Giaquinta and Modica in [8| 



and Weidemaier in [27[, where the local solvability of quasilinear diagonal 
parabolic systems is proved, to show the local existence of strong solution 
to the general transmission problem for isotropic media under less 

restrictive assumptions on the operator B{u) and the parabolicity condition 
of the problem. The main result (local in time existence) is proved by means 
of a fixed point argument based on the Banach contraction principle. 

The paper is organized as follows. In Section [21 we introduce the appropri- 
ate function spaces and recall important embeddings and interpolation-like 
inequalities needed below together with some auxihary results. In Section [3], 
we specify our assumptions on data and structure conditions and introduce 
the precise definition of admissible domains describing the composite body 
under which the main result of the paper is proved. In Section IH we prove the 
existence and uniqueness of the solution to an auxiliary linearized problem 
using the regularity result for elliptic systems in composite-like domains. To 
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make the text more readable, technical details of the proof are collected in 
Appendices [Appendix A[ [Appendix B| and [Appendix C[ The main result 



is proved in Section [5| via the Banach contraction principle. Finally, in Sec- 
tion [6] we present applications of the theory to selected engineering models 
of heat and mass transfer. 



2. Preliminaries 

2.1. Definition of some function spaces and notation 

We denote by W^f = W^'P{ni)^, / > (/ need not to be an integer) 
and I < p < oo, the usual Sobolev space of functions defined in fli and by 
W^~P^^^'^ = W'~^^^'^{dQiy the space of traces of functions from W^'^ on d^e. 
We set = W°'^. Let B be an arbitrary Banach space, then {B)* represents 
its dual. <p'it) indicates the partial derivative with respect to time; we also 
write 

In order to define the concept of strong solution, we will make use of the 
following Banach spaces 

Xe^T := {0; e L-(0,r;L2), e L2(o,r; wf ), 

0"(t)eL2(O,T;L,2), 0(0) = 0} 

and 

respectively, equipped with the norms 

ll</>IU,,T := II0'WIIl-(O,T;L|) + W{t)\\LHO,T;W'f) + II W II L2(0,T;L2) (5) 

and 

Myi,T ■= ll'/^'WllL2(0,T;Lf) + ll¥'(0)||wl-2' (6) 

respectively. 

Throughout the paper, i and m are assumed to always range from 1 
to M and m ^ i, whereas indices i,j = 1,2. Unless specified otherwise, 
we use Einstein's summation convention for indices running from 1 to 2. 
We shall denote by c, ci, C2, . . . generic constants independent on T having 
different values in different places. Let us stress that throughout the paper 
the function C = C{T) depends solely on T and C{T) 0+ for T — i- 0+. 
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2.2. Some embeddings and interpolation like-inequalities 

In the paper we shall use the foUowinK embeddings (recall that f2 is a 
two-dimensional bounded domain) (see (ll. EgI): 



W,^'' ^ Ll < c 11011^1.2 V0 G W]'', 1 < p < oo, 

Wf W,^'^ 11011^1,, <c\\(j)\\^,. \/cf> e Wf , 1< / < 2,p = 2/(2 - /), 
Wj^ < c 11011^. j V0 e Wj^ > 2. 

(7) 

Let us present someproperties of A^t- Assume G A'^t- Using the inter- 
polation inequality [l, Theorem 5.8] 

||0'(t)||L4<c||c/>'(t)||^t||0'(t)||^i' (8) 

we obtain 

< (9) 

For all G Af^^T we have 

||0IU-(O,T;Lf) < C||(/)||^^(0^^.w2,2) < II 0'(t)||^2(o,r;Wf 

< cTV^II^II^^^^. (10) 

Further, combining ([7j) and the interpolation inequality [H, Theorem 5.2] we 
obtain 

U'{t)Ur < cWm^Y^ < c||0'(t)llwv- < c||0'(t)||^l.||0'(t)||i// (11) 

and consequently 

||0'(i)||L3(o,r;L-) < c||0'(t)||i3(o,T;W,i.3) < c||0|U,,r- (12) 

3. Structure conditions and admissible domains 

In this Section, we summarize our assumptions on the problem data and 
specify in detail the geometry of the considered domains. 
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3.1. Structure conditions 

(Al) For every z G M^, B}{s, Z2) and -B|(-2i, s) are increasing functions (with 
respect to s), 5^ : ^ M, such that are bounded on every 

bounded set in for |a| < 3. Further, we denote the matrix 



h y^) - 



(A2) A;^ : X M^^^ — > are continuous and of the semihnear form 

2 

AKr,.) = J]af(r).„ (13) 
1=1 

for all r G and Sj = {s],sf), where si G for z,j = 1,2. Note 
that in f[T^ r stands for u and Sj stands for the vector Vm*. Functions 
a^* : — M, are positive, scalar (due to the assumed isotropy of the 
material) and |(9"a£*(r)| are bounded on every bounded set in for 
I a I < 3. Further we assume 

(14) 

in Q and the ellipticity condition 

in Q with /x^ representing the initial distribution of the unknown fields 

(A3) : — )■ M^, \d"f^{z)\ are bounded on every bounded set in for 
|a| < 2; 

(A4) : r X (0, T) X — )■ is of the form of the Newton-type boundary 
conditions 

gj{x,t,Ui) = a\{u\ - a^{x,t)), 
where are given positive constants and cr : F x (0, T) — )• M^, j = 1, 2, 

a- G 1^2'2(0, T; (Wp^''')*) n W^'\<d, T; Wr^'''). 



/ &]^(/x,)af (m,) + fef (M,)ar(M^) 
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3.2. Admissible domains 

In what follows, we assume that (cf. Figure [1]) 

(i) Q is decomposed into nonoverlapping subdomains Qf, 

(ii) there exists a finite set S C dQ of boundary points such that dfl \S is 
smooth (of class C°°); 

(iii) for every P E S there exists a neighborhood Up and a diffeomorphism 
Dp mapping finUp onto K-pflBp, where /Cp is an angle of size up < n 
with vertex at the origin (shifted into P), 

Kp := {[xi,X2] e M^; < r < oo, < (p < Up} , 

and i?p is a unit circle centered at the origin (r, ip denote the polar 
coordinates in the ( )-plane); 

(iv) the interfaces Tme are smooth (of class C°^), m = 1, . . . , M, m ^ l\ 

(v) there are no cross points of Vmt in ^• 

Let Jv[ be the set of all boundary points A G F = 9fi fl F^^, m = 1, . . . , M, 
i.e. the points where any interface F^^ crosses the exterior boundary d^. 
Further, we assume that 




Figure 1: Admissible domains. 



(vi) for every A G A^, such that A G (9^2 flFmf, there exists a neighborhood 
Ua and a diffeomorphism and DA,m, respectively, mapping QiHUa 
onto /C^_^ n Ba and fi^ fl Ua onto /C^.m H i?^, respectively, where ]Ca/ 
and /Ca,™, respectively, is an angle of size ui and Um, respectively, with 
vertex at A 

^A/ ■= {[xi,X2] e R'^; < r < oo, < (f < ue} 
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and 

^A,m ■= {[Xl,X2] G M^; < r < OO, W£ < V? < Ul + Um] , 

respectively, and Ba is a unit circle centered at the origin (r, denote 
the polar coordinates in the (x'^, a;2)-plane with the origin at A); 

(vii) for every A E M., such that A G dVL fl V^i we have Ui = Um] 

(viii) for every A E Ai, A E dQ fl Tme, we have ua = u;e + ujm = 20;^ < vr. 

Remark. Note that conditions (i)-(viii) incorporate, as a special case, rect- 
angular domain composed of regular rectangles Qe. This serves as a basic 
model for building envelopes, e.g. 

4. Solutions to an Auxiliary Linearized System 

Following the standard methodology of contraction-based proofs, we con- 
sider first an auxiliary linear problem with homogeneous initial condition in 
the form 

in QiT, (16) 

on Ser, (17) 

onr^, x(0,T), (18) 

onr„, x(0,r), (19) 

in n^. (20) 

Assumptions. In (fT6|) -( l20|) are real positive constants, := f3j\x), nf : = 
kI\x) are real positive Lipschitz continuous functions such that 

I5l'l5r4\f>^^l^l±^y inn, (21) 

and the ellipticity condition 

K]^Kf > nfnf in (22) 



i^i h lyjui 



duf 



Uo = u. 



gjix,t) 
j 







U£{x, 0) = 
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Definition 1. Let /, G 3^,,Tand^, G W^^\Q,T- (Wy^'nr)*)nW^i'2(0, T; Wy^';^^) 
Then ui G X^^t is called a strong solution to the system (|T6|) - (l20l) iff 

M ^ i M M 

E / /^^V + E / ^'^^^^^ ■ + E / 

M A/ „ 

= ^ / fiv^dx + J2 / (23) 
holds for every v G W^'^ and almost every t G (0,T). 

Remark. The regularity in time direction naturally imposes two higher or- 
der compatibility conditions on the given functions and in (|T6l) - (fT7ll . 
Namely, the first one requires gf{x,0) to be compatible with (|T7|) while the 
second one roughly says that u'f^(t)\t=o has to belong to appropriate Sobolev 
spaces. This implies additional conditions on f^{x, 0) included in the defini- 
tion of the space 3^^,t- 

Theorem 4.1. Letf, G ye,T, 9e e W''\0,T; (Wy^^;^r)*)niyi'2(o, T; Wy^^;^r) 
and gf{x,0) = on dQ. Then there exists the unique strong solution 
ui G Xg^T to the system f|T6l) - fl20|) and the following estimate holds 



(24) 

To prove Theorem 14.11 we prepare the following definitions and lemmas. 

Definition 2 (Problem (-P/)). Let us define Problem {Pf) by the linear 
transmission system f[TBl)-fl2Ul) with = on dQi fl F x [0,T). 

Lemma 4.2. Let G ye.r- Then there exists the unique strong solution 
Ui G Xi^T of Problem (Pf). Moreover, the following estimate holds 

\\ue\\xt,T ^ 4fe\\yeT- (25) 



Proof. See Appendix A The proof relies on the results for stationary trans- 



mission problem presented in [Appendix C[ □ 



Definition 3 (Problem (Pg))- Let us define Problem {Pg) by the linear 
transmission system flTBl)- fl2(l with = in fi^ x (0,T). 
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Lemma 4.3. Let G W'^''^{<d,T; (W^g'^r)*) n iyi'2(0, T; W^g^r) and the 
compatibility condition gi{x, 0) = on dQ be satisfied. Then there exists the 
unique strong solution U£ G t of Problem (Pg) and the following estimate 
holds 



\U£ 



< c 



lf^llH/2.2(0,T;{Wy„^;^j,)-) 



Proof. See Appendix B Similarly as in the proof of Lemma 



(26) 
we use the 



results for stationary transmission problem presented in Appendix C □ 



Proof of Theorem ^.i, The assertion follows from the superposition principle 
of the solutions to the linear Problems {Pj) and {Pg). □ 



5. Solutions to the Nonlinear Parabolic System 

Definition 4 (Problem (-Po))- Let us define Problem (Pq) by the initial- 
boundary value transmission system with data and structure condi- 
tions satisfying the assumptions (Al)-(A4), see Subsection I3.1[ 

Definition 5. A function u^, such that u'^{t) G L^(0, T; W^'^) and u'f it) G 
L^(0, T; L^), is called a strong solution of Problem (Pq) on (0, T) with initial 
data G Wf'^ iff 

M „ M 

^^-^ Jdnenr Jui 

holds for every v G W^'^ and almost every t G (0, T) and 

Ui{x,0) = fj,£{x) in rii. 

Theorem 5.1 (Main result). Let the assumptions (A1)-(A4) be satisfied. 
For a given ^t£ G W^'^, which is supposed to be compatible with ([2]) -([3]), 
there exists T* G (0, T] and a function such that U£ is the strong solution 
of Problem (Pq) on (0,T*). 

Proof of the main result is postponed to the end of this section. We start 
from a related problem with homogeneous initial condition. To that end, let 
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Ui be the strong solution of Problem (Pq) on (0, T), — fj,^ + y^. Then 
& Xj^^T and the following equation 



M ^ i M 



+ E / + y'i - d5 = E / + y^) d^ 



£—1 J dCl^nT JQ( 

rl,2 



holds for every v G W ' and almost every t e (0,7"). This amounts to 
solving the problem with shifted data 

br{x,yc) = bPiy^ + He), 

ft{x,y^) = f^ini + y^). 

We often omit the argument "a;" writing shortly af'{y^) instead of a;^*(a:, y^), 
^/{ye) instead of ^^/{x^y^) and f^{yi) instead of /^{x^y^). 

Definition 6. Define the operator J^f : X£^t ye,T given by 

E/ J^M-vdx ^ f^f (brm-bri<t>e))^v^dx 

M 



+ E / («f (0) - «f l'^^)) • V^^- da; 

M 

- E / «i'(</>^) V/^^ • V^;^' da; 

M „ 

+ E / fii^e)^'^^^ (27) 



which holds for every v G W^'^ and almost every t G (0,T). 

Remark. Let W£ = fi^ + yf . The function ti£ is the strong solution of Problem 
(Pq) on (0)^) with initial data fj,^ G W^'^ iff for y^ G the following 
equation 



E / ?(0)i^^^da; + E / af{0)Vyi-^vUx + J2 / d5 
+ E/ al{^il-a^)v^dS = J2 ^{ye)-vdx 
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holds for every v G W^'^ and almost every t G (0,T). 

Before proceeding to the proof of the main result of the paper, we prepare 
some auxiliary lemmas and propositions. 

For a fixed R> define the closed ball Br{T) C X^^t 

Br{T) :={0gA',,t; \\ct>U^r<R]- 



Lemma 5.2. Let (p^ G Br{T). Then 

W^MWy,, < CiC{T) + \\'t>i\\XT + II'^^IU^.T ) + C2, 



(2^ 



where the function C{T) — j- 0+ for T — 0+ and the constants ci, C2 > 0, both 
independent of 4>^, do not depend on T. 

Proof. The proof is rather technical. To derive the estimate (128|) we exten- 
sively use the embeddings and estimates (17I)- (fT2|) . First, for all G Br{T) 
we have 



11^(0. 



Wily,. 



< 



S,(O)-6,(0,))<^;(t) +||V-[(af(O)-af(0,))V0^] 



yi 



yi,T 



+ ||V- [af(0,)V/i^] 



\yij 



ye,T 



(29) 



Now, we have to estimate each term on the right-hand side of (l29l) . Succes- 
sively, we use ([9]) and (Al) (see Subsection 13. ip to estimate the first term: 



< 



yt- 



S,(o)-6,(0,)U;'(t) 



+ 



m$^Mm\)\t) 



< Cl||(/),||Lcx,(Q^^)2||(/)"(t)||i2(Q^^)2 

IT) 

< CiTV2 110^11^^^^ +c2TV4||0;(t)||2^^^^^^^^^ 



il/4 



(30) 



Similarly, estimating the second term in fl29|) in the norm of the space yi^T 
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we arrive at 



|V-[(af(O)-af(0,))V0a||j;,^< 



+ 



+ 



+ 



VU[)'{t)]-Vct>\ 



dx^d(l)\ 



+ II V[af (0) - a\\x, 0,)] ■ V(0y'||^,(Q^^), . 

(31) 

The first integral on the right hand side in f l3ip can be estimated 



< C 



|0Ki)l'|V0/da; dt 



\Jni 

T 







< c||0^(t)||^2(o,T;L|) II<^^IIl°°(0,T;W,''*) 



and applying ([9]) we obtain 



(32) 



L^QtT? 



Similarly 



^S^V[(4)'(t)] ■ 



< c 



\Vct>[{t)\^\Vct>,\'dx ) dt 



VJOf 

T 



^IIl°°(o, 





< c||(/)^(t)||^2(o,T;W,i'3) II<^^IIl-(0,T;W,^''') 

< cT^/3||0;(^)||5.,.^w,3^^l^2"" 
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and using f lT^ we get 



(33) 



Similarly, the third term in (13T]) can be estimated as 



; We) v^) 



dx^d(j)\ 



(34) 



Further 



< c 1^ i0Kt)nA(/>/ dxj dt 

< C / \\(f)[it)\\l^ ||0J^2,2dt 

Jo 

< c||0^(t)||^2(o,T;Lf ) II^€IIl-'(0,T;W^'2) 

< cTV3||0;(t)||^3 tV2||0^||^^^^.(35) 



Now dSSD and imply 



90 



TOj) yields the estimate 



^^^('^^\4)'(t)A0^ 



(36) 



< cT'l^4>e\\l,r. (37) 

where we have used the Lipschitz continuity of a^. In the similar way one 
can deduce 

||V[af(O)-af(0,)]-V(0^)'(t)||^,(g^^), 

< Ci||V(/)£(t)||2,2(o,r;L|) + C2 II0^IIloo(o,T;WJ''^) L2(o,T;WJ'^) 
< CiTV6 ||0;(t)||^3(o,T;W-) + II0K^)IIl3(O,T;W-) 

<ciTi/l0,|U,,, + C2r^/il0,||^^^^. (38) 
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Finally, the estimates fl^ -f l55]) imply 

\\V-[{af{0)-a\\4>,))Vct>\]\\y^^^ 

< cc{T) + \m%^^ + ii0,iu,,) , (39) 

where the function C(T) — 0+ for T — > 0+ and c is independent of T. 

Further, estimating the third term on the right hand side in (129|) one 
obtains 



< 



+ 



+ 



LHQit) 



+ ||V- [af(0)V/xa||^,.. (40) 
Estimating each term on the right hand side we arrive at 



■(0^)'(t)V0^V/i^ 



<c \\m\\lf\\'^Mh\\Vfi,\\isdt 

/n £ £ £ 

< cT3/1l0;(t)||is(o,T;L|)ll'^.lli.(o,T;W^«) HVM.IIlB 

<cT'/'\m%^r, (41) 



where we have used the estimate ||</)£|P^ l^8^ < T^^'^W'PeW'xt ■ Further 



< cj^ \\'^m)\\^\\'^^^ir^dt 

< cTV3||0;(t)||^3 3^||VMj?.a 



(42) 
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< c 



\(p'eit)\\lr \\^^^e\\h 



< cTl/i0;(t)||i3(o,T;L?=)l|A/^£llL- 



(43) 



and finally 



< cTV3||0;(t)||23^ ||V;.,||2 



(44) 

The inequalities (HUll - lH^ yield the estimate 

II V ■ [af (0,)V/i^] 113,^^^ < CiT^/^ 110,11^^^^ + + (45) 

where the constant C3 in (145|) bounds the last term in (140|) . Finally, taking 
into account (A3), the source term can be estimated as 



yi,T 



< cT^/^||0;||i3(o,T;L|) +C2 



(46) 



Altogether, the estimates ([30]), ([39]), (jlS]) and (jlH]) yield the inequality (|28]). 

□ 

Lemma 5.3. There exists a nondecreasing function c{R) (c{R) does not 
depend on T , (f)^ and (f>^) such that for all 4'g,4'g G Bfi{T) 



||jr(0,) - Jr(0,)||3;,,, < C(i?) C(T)||0, - 0,|U,_„ 
where the function C{T) — ^ 0+ for T — )■ 0+. 



(47) 



Sketch of the proof. Similarly to Lemma 15.21 the proof is rather technical. 
Therefore, we only sketch the procedure and omit the detailed derivations. 
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First, we estimate 

< 11(6,(0) -6,(0,)) cf>',{t) - (6,(0) -6,(0,)) 4>'e{t) 
V- [(af (0) - af (0,))V0^] - V- [(af (0) - af (0,)) 
V • [af (0,)V/.^] - V ■ \af{4>e)Vfi: 



ye- 



+ 



+ 



ye: 



ye.T 



ye: 



(4J 



The right hand side in fBHl) can be further estimated by 



6,(0,)- 6,(0,) U;(t) 



+ 

ye.T 

6,(0,) (0;(t)-0;(t) 



6,(0) (<t>'eit) - 0,(t) 
+ 

V(0)V 
V 



ye,T 



V- 

+ 



ye: 



+ 



+ 



V- 



af (0,) - af (0,)) V/i; 



af(O)V(^0^-^ 
af(0)-af(0,))v0; 



ye,T 
ye.T 



(49) 



Estimating each term in using the same arguments as in the proof of 
Lemma [5.21 and the assumptions (Al)-(A3), see Section [XT| one obtains the 
inequahty 

II jr(0,) - ^(0,)b,, <c,{R' + R + 1) C(T)||0, - 0,|U,, (50) 

V ' 

c{R) 

for all 0,, 0, G Sij(T). Now ([50]) yields dH]). □ 
Using Definition |6] and Theorem 14.11 we can formulate the following 

Definition 7. Let =Sf : Xe^T '^e,T be an operator such that =^(0,) = 
if and only if 



M 



e=i 



dy\ 



M 



M 



a". 



1=1 

M „ 

£=1 



holds for every v G W^'^ and almost every t G (0, T) and ?/,(a;, 0) = in Qg. 
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Proof of Theorem \5.1\ (Main result). The proof of the main result is based 
on the Banach fixed point theorem. Lemma 15.21 and the estimate fl2^ imply 
the inequality 

< C2C{T) (110,11^,^ + 110,11^,,^ + ||0,|U,,,) + K2 (51) 

for all 0^ G Bfi{T), where Ki and K2 are positive nondecreasing functions 
with respect to T and independent of and the constants ci,C2 are inde- 
pendent of <pi and T. Further, linearity of (fT6l) - (l20|) . the estimate ( HM and 
Lemma 15.31 imply 

< c{R) C(T)||0, - for all 0„ 0, G i?^(T), (52) 

where c{R) is some nondecreasing function and C(T) — )• 0+ for T — i- 0+. 
Now 0511) and 0521) imply that for sufficiently small T* G (0, T] there exists 
R > such that =Sf : A^^* — ^ '^^.t* niaps Br{T*) into itself and =Sf is a 
strict contraction in Bfj{T*). Hence, using the contraction mapping principle 
we have the existence of a fixed point G B^{T*) C X^^t*-, such that 
■^{.Vi) = Vi- Vi is uniquely determined in the ball BuiT*). Set = fx^ + y^. 
By Remark O the function is the strong solution of Problem (Pq) on 
(0,T*). □ 



6. Applications 

In this Section, we present examples of the coefficients of the parabolic 
system ([T]) related to models of transport in porous media. Note that for 
brevity, we omit the subscript i and the dependence of all variables on x and 
t in what follows. 

All available engineering models of simultaneous heat and moisture trans- 
fer possess a common structure, derived from two balance equations of heat 
and mass joj: 

dH . ^ dM 

^ = -V-jQ + g, -^ = -V-j^, (53) 

where H (Jm~^) is the specific enthalpy, M (kgm~^) denotes the partial 
moisture density, Q (Jm~^s~^) stands for the intensity of internal heat sources 
and Jq (Jm~^s~^) and (kgm~^s~^) are the heat and moisture fiuxes, 
respectively. This structure is also refiected in the choice of the unknowns u, 
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which consist of the temperature v} = 9 (K) and a quantity related to the 
moisture content. 

Individual models are then generated by the choice of the second state 
variable u"^ and of the individual components in system (l53i) . In Sections 16.11 



and 16. 2[ following the expositions of Dalfk et al. and Cerny and Rov- 
nanikova j3| , we briefly introduce two such representatives due to Kiessl 12 



and Kiinzel [18|, simplified by assuming that freezing of water in pores has 
a negligible effect. An interested reader is referred to 0, 0, [3] for addi- 
tional discussion of the models and to ^ for details on the terminology used 
hereafter. 

6.1. The Kiessl model 

The enthalpy term in the Kiessl model is postulated in the form 

H = pqCqO + pu,Cww6, (54) 

where po (kgm~^) and cq (Jkg~^K~^) denote the partial density and the 
specific heat capacity of the dry porous matrix, and Cyj are analogous 
quantities for water and w (-) is the relative moisture content by volume. 
The heat flux follows from the Fourier law for isotropic materials 

jQ = -\{w,e)we, (55) 

with A (Jm^^s^^K"^) being the state-dependent coefficient of thermal con- 
ductivity. 

The moisture balance is based on the moisture density provided by 

M = pyjW + (e - w)ippp^s{0) (56) 

where e > w (-) denotes the porosity, ip is the relative humidity and pp^s < 
Pw (kgm~^) is the material-independent partial density of the saturated vapor 
phase, given as a smooth increasing function of 6. Assuming again isotropy 
of the material, the corresponding moisture flux is expressed in the form 

= - {D^iw, e)Ww + D^{w, e)Vip + Deiw, 6)^6) , (57) 

where Dyj (kgm^^s^^), (kgm~^s~^) and Dq (kgm~^K~^s~^) denote material- 
specific diffusion coefficients, which need to be determined experimentally. 
Finally, the internal heat sources 

Q = L, {D^{w, e)Vw + D^{w, e)Vip - ^ [(e - w)ippp,si0)]) (58) 



19 



quantify the influence of phase changes in pores by means of the latent heat 
of evaporation of water L„ (Jkg~^). 

To close the model, Kiessl in [l2| related the auxiliary variables w and 
ip to a dimensionless moisture potential $ = via monotone, material- 
dependent, transformations 

w = /($), ip = 9{^), (59) 

satisfying /(O) = g{0) = and ^(0) = 1. In particular, / denotes the 
sorption isotherm, whereas g reflects the pore size distribution. By employing 
these identifies, the individual coefficients in ([1]) receive the form (cf. j^]) 

= poCoe + p^c^gi^)e + L,ie- fi^))gi^)pp,,i9), (60a) 
= p^/($) + (e-/($))^($)pp,,(e), (60b) 
a'' = A(/($),^), (60c) 

= L^(^DUm),0)^^ + DMi^),0)^y (60d) 

= De{fi^),e), (60e) 

DUm), 0)^^ + D,{f{^), 6)^. (60f) 



6.2. The Kiinzel model 

In the Kiinzel framework, the heat balance is described using identical 
expressions for the enthalpy fl5^ and the heat flux ( 155|) as previously. In 
addition, the moisture density is simplified into 

M = p^w (61) 

and the moisture fiux attains a form 

j„ = -{D^(^,e)\/cp+^-^v i^psm), (62) 

p 

in which (kgm^^s^^) stands for the liquid conduction coefficient, 6 (kgm^^s^^Pa^^) 
is the vapor diffusion coefficient in air, p (-) is the vapor diffusion resistance 
factor of a porous material and Ps (Pa) is the vapor saturation pressure. This 
yields the heat source term given by 

Q = L^V ■ (^^V i^Psmy (63) 

The relative humidity is chosen as the second unknown, u"^ = cp, and is 
used to express the associated volumetric moisture content in the form 

w = h{p), (64) 
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where h is a monotone moisture storage function with h{0) = 0. Altogether, 
the coefficients in system ([T]) read as 



= pqCqO + p^c^,h{ip)9, (65a) 

= p^h{ip), (65b) 

a" = X{h{^),e) + Lj-^^^, (65c) 

p dtJ 

= lJ-^pM, (65d) 



p 

p ^ d9 



a = ip , (65ej 



a^^ = D^{^,9) + ^-^ps{9). (65f) 

5. 5. Structure conditions (Al) and (A 2) 

The structure conditions (Al)-(A2) closely reflect the physical constraints 
on the underlying transport models and experimental observations. Con- 
cretely, the model parameters (such as e.g. /($), g{^) and Pp,s{9) in the 
Kiessl model, or h{(p) and Ps{d) in the Kiinzel model) are obtained by fltting 
smooth functions to experimental data, determined for a limited range of 
state variables. The required regularity and boundedness of coefficients B, 
a^^ and positivity of a^'' is therefore ensured. The increasing character of B 
is consistent with the fact that both the speciflc enthalpy H and the moisture 
density M increase with an increasing temperature and the moisture-related 
variable, respectively. The ellipticity condition (IT^ is satisfled due to the 
fact that the Soret- and Dufour-type fluxes, quantifled by a^^ and a^^, are 
dominated by the diagonal contributions a^^ and a^^, see also i. Therefore, 
any physically correct form of a-'^ must meet this condition. The validity of 
the assumption f[T4|) then follows from the same physical reasoning. 



Appendix A. Proof of Lemma 14.21 

Discretize fl25]) in time and replace u'^{tn) by the backward difference quo- 
tient d^^{w£)n = [{wi)n — {wi)n-i]/h, where /i > is a time step. Suppose 
r = T/h is an integer. For simplicity, let us write Wi = {wi)n, fi = {fe)n- 
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We have to solve, successively for n = 1, ■ ■ ■ , r, the steady problems 

M M 

+ J2 d5 = ^ / fi dx (A.l) 

which hold for every v G W^'^ and {w£)q = in Qi. Test (lA.ip by [f^,f^] = 
[nj^ip^, i^Y'^'^] ^^"^ define the bilinear form 2t(iu, 

^ Af „ Af „ 

2t(w;,, ^) = -^J2 da; + ^ / /tf /tf V^i ■ Vy^^' do; 

Af „ 

+ X] / w^^^d^, p=l,2, p^j, 

for every (p G W^'^. Set = W£. Now, fl?T]) . ([22]) and the Friedrichs 
inequality yield the W^'^-ellipticity, i.e. there exists c > such that 

c\\w e\\l,^,2 <\Qliwe,We)\ (A.2) 



for all We E W^' . Using the Holder inequality and the standard trace 
theorem one obtains the continuity of 21, i.e. the inequality 

\^{we,ze)\ < c\\we\\^i,2\\ze\\^i,2 (A.3) 

1 2 m 

which holds for all we,ze G W^' and for some positive constant c. The 
linearity of 21 : W^'^ (W]'^)* is obvious. Hence, for every G C 
(W^''^)* the Lax-Mil gram theorem yields the existence of the weak solution 
Wf G W]'^. To get higher regularity results (with respect to time), define 
Wi G W^^'^ by (jXU) and test (jM]) by [v\v^] = \Kfdt^w\, nfd^^wj] to get 

A/ 

Af „ A/ „ 

J fie JdQef^r 
M „ 

= E / P = 1, 2, p ^ J. (A.4) 
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Denote by 



+ Kl\fV{wl)n ■ V{w^^)n, n = 1, . . . , r. (A.5) 

Now we can estimate 

(!>'^[{Vwe)n] ■ {{Vwe)n - (Viu^)n-i) > $4(Viu^)„] - ^e[{Vwe)n-i] 
because is convex. Thus, using Young's inequality, one obtains 

i=l i=l •'^i 



1 



£=1 



PI 1 



[w 



^ - iiK)n-ir 



M 



M 



<C{e)f2\\fi\\U + ^c,f2 f Wd^'weWl.dx, p=l,2,p^j, (A.6) 
e=i e=i 

with some arbitrarily small constant e. Note that (12T]) yields the estimate of 
the parabolic term 

M „ M „ 

where c depends on f3j^ and k-^*. Based on the estimate (lA.6p . the same way 



as m 



22I . Proof of Theorem 8.16] we can prove the existence of the solution 
Ui G L°°(0,r; W^^'^) L°°(0,T; L^) with u'g{t) G 1.2(0, T; L^) and the 



estimate 



I'"K^)IU2(0,T;L2) + II'"£||l°°(0,T;L2) < c|| || 2,2(o,T;L2) ■ 



(A.7) 



Now we can proceed as in 10||. Rewrite the system (|T6i) - (l20l) in the form 



-V • (nfVul) 


- 

— ^ I 




= 


K 




i drii '^rn dtlm 


= 



a dui 



in QiT, 

on SiT-, 

on Tmi X (0,T), 

on Vmi X (0,T). 



(A.8) 



Since u'{t)i G 1.2(0, T; Uf) we have F^g{x,t) G for a.e. t G (0,T). Accord- 
ing to results for stationary transmission problem (see Appendix C Corol- 
lary [1]) we have Ui{t) G W^'^ and the estimate 



Ue{t)\\^22 < c||/^(t)i|L2 



(A.9) 
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holds for almost every t G (0,T), where the constant c is independent of t. 
Raising f lA.Qp and integrating both sides in the inequality (lA.Op with respect 
to time, we get Ui G L^(0,T;W^'^) and taking into account the estimate 
(1A.7|) we arrive at 



\\f^'ei't)\\L^O,T;Vj) + ll'"dlL2(0,T;W2-2) + II ""^1 (0,T;L2) < c|| || £2 (o,T;L2) • (A. 10) 

Now taking the derivative of (fT6|) - (!20|) with respect to time, considering 
fi ^ yi,T (including the compatibility condition on f^{x,0)), we conclude 
u'lit) G L2(0,T;L2), u'g{t) G L^{0,T;Wf^) D L'^{0,T;Lj) and the estimate 
fl25|) follows. The linearity of Problem (Pj) and the estimate fl25|l yield the 
uniqueness. 

Appendix B. Proof of Lemma 14.31 

Let G W^'^{0,T; (Wg(J'^r)*). We are looking for the solution of the 
problem defined via 

M „ M 

+ E/ ^■e>^d^ = 5^(G^'(t);i;)( 1,2 ,2 (B.l) 

to be satisfied for every v G W-*^'^, almost every t G (0,T) and ^^(a;,0) = 
in Qi. The duality of {Gi, v)iy^i,2y ^1,2 corresponds to 

M „ 

{Ge; t;)(w,^. wM = J2 9l dS. (B.2) 

Approximate fIB.ll) in time by discretization and replace ^^(tn) by the back- 
ward difference quotient d^^{wi)n = [{wi)n — {we)n-i\/h, where > is 
a time step. Suppose r = T//i is an integer. Let us write Wi = {wi)n and 
test fIB.ip by [t>^,t>^] = [nYip^, nj'^Lp'^]. We have to solve, successively for 
n = 1, ■ ■ ■ ,r, the steady problems 



M ^ M 

PJ J J 7 



M 



E/s:f((G"'(t„))^,;^^')(^i,2).,^i,2, p=l,2, p^j, (B.3) 
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to hold for every <^ G W^'^ and {10^)0 = 0. Test (IB.Sp by = {wi)n to get 
the estimate 



Cl 



M , 



Wi)m\ dS" 



1/ n 



1,2, p=l,2, p^j. (B.4) 



= 1 771=1 



Now we can proceed as in 22, Proof of Lemma 8.6, Proof of Theorem 8.9] 
to prove the existence of the weak solution G L^(0,T; W]'^) (as the limit 
of Rothe sequences) with the estimate 

Let us note that stands for u'^{t). Hence u'p{t) G L2(0,r;W^'^) 
L2(0,T;L2). Further, let G ^^^'^(O, T; Wgg^r). The standard theory 
for parabolic problems yields u'i{t) G L°°(0, T; L^). Using the same proce- 
dure as in Appendix A and according to results for stationary transmission 
problem (see Appendix C[ Corollary[T]) we conclude u'pit) G i^^(0, T; W^'^)n 
L°°(0,T;L2) and (combining with flB31) ) 

ll'""WllL2(0,r;L2) + II'"£(^)IIl2(0,T;W2'") + II ^K^) II L°°(0,T;L2) 

< c (||»,|li^2,2(o,T;(wy„^;^,).) + ll»JH^i.2(o,T;wy^^;^,)) ' (B-6) 
The linearity of Problem {Pg) and the estimate (IB.SP yield the uniqueness. 
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Appendix C. Transmission problem for elliptic systems in a multi- 
layer structure 



The boundary transmission problem for the elliptic system for the sub- 
domain is formulated, in the expanded form, as 



-V 
-V 



-11/ 



.21, 



X 



X 



duj 



+ e 
+ e 



12, 



22, 



X 



X 



-111 



-21, 



duj 
duj 



X 



X 



duj 
du] 

arte 



e 



V- 
V- 

12 



{ef{x)Vui) 
{ef{x)Vui] 



+ e 



22, 



X 



X 



duj 
duj 

arte 



+ e 
+ e 



X 



X 



8"m 

an„ 
an„ 



+ e 
+ e 



- dui, 



X 



X 



drirr. 



f! 






m 
in 

on 

on 
on 
on 

on 

on 



mlt 



nr, 
nr, 



(C.l) 

Here we assume that the problem (IC.ip is elliptic and has a unique weak 



Further we consider 



solution Ui G W]'^ for G and G W](?'^. 
that s-'l{x) are positive Lipschitz continuous functions and a\ are prescribed 
constants. 

Elliptic boundary value problems in cornered plane domains are exten- 



sively investigated in the literature, see e.g. p, [13 



14, 15, 171. The behavior 



of local solutions of g eneral linear and semilinear transmission problems is 



studied in 



23|,l2J, 
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We adapt the general framework stated in the liter- 
ature to calculate the regularity of the plane transmission problem for the 
elliptic system of equations fIC.ip . 

It is known (cf. [l7|, l25|) that, in general, the boundary singularities may 
occur near corner points at the boundary dVL, the points at the boundary 
where the boundary conditions change their type, the crossing points of inter- 
faces, corner points of inclusions or points, where the interfaces Tmi intersect 
the exterior boundary of the domain f2. Taking into account the assumptions 
on admissible domains introduced in Subsection 13.21 only the points where 
the interfaces T^i intersect the exterior boundary are of importance in our 
analysis of vertex singularities. Hence, let M. be the set of all boundary 
points y4 G r n Tmii m, i = 1, . . . , M. As well known, the local regularity is 
valid outside an arbitrarily small neighborhood of the points A E M.. Hence 
it suffices to prove the regularity for the solution ue with small supports. For 
solutions with arbitrary support the assertion then can be easily proved by 
means of a partition of unity on Q. Let A be an arbitrary point from the 
set Ai and let the support of U£ be contained in a sufficiently small neigh- 
borhood U{A) of the point A. Let D be a diffeomorphic mapping fl r\U{A) 
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onto }Ca(^Ca, where JCa is an angle with vertex at the origin (shifted into A) 
and Ca is a unit circle centered at the origin. Considering a zero-extension 
outside Ca, a new boundary value problem with "frozen coefficients" ef{0) 
is defined in an infinite angle Ka = /C^ U /C^+i and the problem of regularity 
of the solution coincides with the original problem near the corner point A. 
One applies certain regularity theorem for the boundary value problem in 
an infinite angle JCa- Hence, following 23|, Chapter 3 and 4], we localize the 
boundary value problem fIC.ll) (see Fig. IC.2I) . i.e. identify the origin of 
coordinates with the point A, "freeze" the coefficients ef{0) and multiply 
the weak solution by a cut-off function ?7(|ic|) G C°°(R^), < ''^(la;!), given 
by 

1 for \x\ < e/2, 
for |a;| > e. 



ri{\x\) 



Denote Wr 



m 



+ 1, and consider the following problem with 



dICe n dlCi 




dICi+ 



Figure C.2: Localization principle. 

frozen coefficients (now written in compact form using the Einstein summa- 
tion convention for indices i and j running from 1 to 2) 

-£^^(0)A< = Fi in /C„, 

- = GL on 9/C™nr, 



vol 



8w 



wl^^ on S/C^nS/C^+i, 



dUp 







on S/CfflS/C^+i, 



where 



pi 

m 



Q3 



dxh dxi 



9rnV ^rn'^rnV- 



(C.2) 
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The regularity of in a neighborhood of the point A is determined by the 
smoothness of near A. Using polar coordinates (r, u) in flC.2p we arrive 
at the system 



' -^^^(0) + + ^_^) = Fl in S^^m = lA + 1, 

5pO)g(r,0) = G^(r,0), 



where and 5*^+1, respectively, is an infinite half-strip 



(C.3) 



Si = {{r,uj) : r eM.+ , < u < Ue} , 
Si+i = {(r, w) : r G IR+, Ue < u < w^+i} , 

respectively, and Wm{r,u) = Wm{xi,X2), Fm{r,u) = F„(xi,X2) and G^(r,w) 
Gm{xi,X2), m = + Substituting r = e^, we get the system of equations 

-^ml^ + 'S') = F^inS^^, m = £,£ + l, 
^f(O)Sf(e,0) = GJ{^,0), 
4i(0)^(e,c^m) = g+i(e,^m), (^-^^ 

where S(, and S'^+i, respectively, denotes an infinite strip 

Si = {{^,u):^eR,0<uj<uji}, 
Si+i = {{^,u) : ^ eR, Ui <uj < Ui+i} , 

respectively, and w{C,, u) = w{xi, X2), Fm{C., u)e~^^ = Fm{r, u), Gm{^, uj)e'^ 
Gm{r, u). We apply the complex Fourier transform J^^^\ with respect to real 
variable ^ G M, 

[J-5^A0(O](A) = 0(A) = ^ / 0(Oe"''^ de, A G C, 
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Si 















m 






'///////////////////// 


J. 








Figure C.3: The infinite half-strip S* = S'^ U S'c+i and the infinite strip S = SgV} S'^+i- 



to transform flC.4p to the one-dimensional problem 

92 



5f(0) (iA)2u;i + 



(A,0) 



F/ for uj G (0, 

for UJ e {ui,uje+i] 



(C.5) 



sf{0f-§iX,ue) 



-e+i 



:of-%^{X,ue) 



with complex parameter A. Solvability of the parameter dependent boundary 
value problems were studied in [17|. Roughly speaking, the operator pencil 
2t^(A), associated with the parameter dependent boundary value problem 
fIC.Sp . is an isomorphism for all complex parameters A G C except at certain 
isolated points - the eigenvalues of 21^ (A) (for precise definition of eigenval- 
ues and corresponding eigensolutions we refer to monograph 17|). As well- 
known, the regularity of the weak solution U£ G W^'^ (the existence of the 
strong solution, i.e. whether or not Ue G W^'^) depends on the distribution 
of the eigenvalues A of the operator pencil ^a{X) in the strip Im A G (—1, 0). 
This assertion is expressed by the following theorem, which is a classical 
result, see 0, 0, Q, Q: 



Theorem Appendix C.l (Regularity theorem in an infinite angle). Let 
We G W^'^(/C^) be the uniquely determined weak solution of flU.2p . G 
hjlJCe), Ge G W]{?'^(/C^). // the strip ImX G (—1,0) is free of eigenvalues 

t2,2/ 



of the operator pencil 'QIa{X), then We G W^' (JCe) and 

(\\Fihj(K:,) + \\Gi 



1/2,2 



Proof. Theorem Appendix C.l is a consequence of [171 . §7, Theorem 7.5] and 
the expansion of the solution (itI. §7, (7.10)]. □ 



29 



The characteristic determinants and the distribution of the eigenvalues of 

Every Aq G C such that ker 21a(Ao) 7^ {0} is said to be an eigenvalue of 
21^ (A). The distribution of the eigenvalues of the operator 2lyi(A) plays crucial 



role in the regularity results of the solution, see Theorem Appendix C.l In 
order to calculate the eigenvalues of the operator pencil 21^ (A) we look those 
A, for which there exists the nontrivial solution of the system f lC.Sp with the 



vanishing right-hand side. The general solution [e^,e^] of the homogeneous 
equations 

-4^(0)((iA)^^^ + ^)-.^^(0)((iA)^e^ + ^) = 
-^^^(0)((iAre^ + |§^)-e^^(0)((iA)^e^ + ^) = 



m 



+ 1, has the form (recall that e^^ is a positive definite matrix) 

(C.6) 



ej = Ci cos(iAa;) + C2 sin(iAa;) for u e (0, ug), 

e| = Cscos^iXu) + C4sin(iAco') for u G (0,c<;^), 

= C5 cos(iA(X') + C6sin(iAu;) for u G {ue,coe+i), 

e|+i = C7Cos(iAci;) + C8sin(iAco') for u G {ui,ue+i). 



The eigenvalues of 21^ (A) are zeros of the determinant Da_{X) of corre- 
sponding matrix of coefficients Ci, . . . , Cg (substituting the general solution 
(IC.6p to the corresponding boundary conditions and transmission conditions, 
respectively, we get the homogeneous linear system of eight equations with 
unknowns Ci, . . . ,Cs). Computation of -Da(A) leads to the transcendent 



equation 
where 



D^{\) = D]li\)Dfi\) - D\'i\)Dl\\) = 0, (C.7) 

-^a(A) = £^"^(0) sin(iAa;^) cos[iA(a;^+i — w^)] 

+£^+1(0) cos(iAu;£) sin[iA(w^+i - ue)], 
= 5^^(0) sin(iAci;^) cos[iA(a;^+i — cj^)] 

+4+1(0) cos(iAa;£) sin[iA(w^+i - cu^)], 
-^a(A) = £^"'^(0) sin(iA(X'^) cos[iA((X'^+i — cj^)] 

+£^+1(0) cos(iAa;£) sin[iA(w<?+i - ue)], 
-^a(A) = £^(0) sin(iAa;^) cos[iA(a;^+i — w^)] 

+£^+1(0) cos(iAco'£) sin[iA((X'^+i — tug)]. 
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The roots of the equation Da{X) = are the eigenvalues of 21^ (A). 

Taking into account the special type of geometry, namely ue = Ui+i/2, 
(1C.7I) simplifies into 

-(.f (0) + e]lM)ief{0) + ejUO))] sin(2iAu;,) = 0. (C.8) 

Since both matrices, e^^ and e^^_^_-^^, are considered to be positive definite, we 
get 

sm{2iXuJe) = 0, (C.9) 

from whence we obtain 

iA = , keZ. 

2uji 

Now it is clear that for Ui G (0,7r/2] there are no roots of the equation 
Da{X) = such that ImA G (-1,0). 

1 2 

Corollary 1. Let ui G W^' be the uniquely determined weak solution of 
flUJD . G Lj, G Wj(?'l Since the strip ImA G (-1,0) is free of 
eigenvalues of the operator 21^ (A), we have Uf G W^'^ and 

||*tdlv\A2,2 < c (^II/^IIli + 1 1 c/£ 11^1/2,2) . 
Sketch of the proof. The assertion follows from Theorem Appendix C.l and 



the determinant equation flC.7l) . flC.7l) implies that for w^+i < n, ue = a;^+i/2 
there are no eigenvalues of the operator 2l^(A) in the strip 1mA G (—1,0). 
Hence W£ G W^'^(/C£). Now consider the weak solution U£ G W]'^ of (IC.ip . 
Let A^^ be the set of all boundary corner points A G dfli, i = 1, . . . , M. We 
have 



= I i - ^ ?7a I + ^ rjAUi (C.IO) 

V AeMe / AeMe 

The regularity of the first term on the right-hand side follows from the interior 
regularity. The smoothness of the second term follows from the regularity 



result in an infinite angle ICi, Theorem [Appendix C.l[ □ 
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